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FOREWORD 

Tbls document presents the results of work performed 
by the Advanced Fluid Mechanics Section of the Lockheed 
Missiles & Space Company, Inc., Huntsville Research & 
Engineering Center for NASA-Johnson Space Center, Houston, 
Texas. 

The contents of this document represent partial ful- 
fillment of the requirements of Contract NAS9- 16256 with 
NASA-Johnson Space Center. The Contracting Officer’s 
Representative for this study was Mr. Barney B. Roberts, 
EX43. 
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1. INTRODUCTION AND SUMMARY 

The overall objective of this study was to deliver a nozzle -plume flow- 
field code that has capabilities which do not presently exist in a single com- 
puter code. The RAMP code (Refs. 1 and 2) which was developed by Lockheed 
under government funding was chosen as the basic code from which to work. 
The basic RAMP code is of modular construction and had the following capa- 
bilities: 

• Two-phase with two-phase transonic solution 

• Two-phase, reacting gas (chemical equilibrium reaction 
kinetics), supersonic inviscid nozzle/ plume solution 

• Operational for inviscid solutions at both high and low 
altitudes. 

During the course of the study the following capabilities were added to the 
code: 

« Direct interface with JANNAF SPF code (Ref. 3) 

• Shock capturing finite difference numerical operator 

o Two-phase, equiiibrium/frozen, boundary layer 
analysis 

e Variable oxidizer-to-fuel ratio transonic solution 

• Improved two- phase transonic solution, and 

• Two- phase real gas semi-empirical nozzle 
boundary layer expansion. 

Most of the above capabilities already existed in other computer codes. 

These codes were incorporated into the RAMP code to enhance its use- 
fulness. 
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This report , discusses the capabilities which have been added along 
with how to use these capabilities. Detailed descriptions of the theories 
are left to the references from which the codes were taken. 

Also included in this report is a detailed input guide for the RAMP 

code. 
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2. REACTING AND MULTIPHASE (RAMP) COMPUTER PROGRAM 

A precise knowledge of local flow properties in nozzles and exhaust 
plumes is necessary for performance, radiation, attenuation, heat transfer 
and impingement analyses. The reacting and multiphase (RAMP) computer 
program is designed to give detailed flowfield information in the supersonic 
region of a reacting multiphase two-dimensional or axisymmetric flow field. 
The boundaries of the flow field may be solid such as in a nozzle or "free" 
such as in a plume. The analysis may be utilized therefore to predict per- 
formances as well as plume characteristics of a given engine system. A 
printed record of the program results is given for user inspection while a 
binary tape is provided for subsequent manipulation by other analyses. 

The flow of a gas -particle mixture is described by the equations for 
conservation of mass, conservation of momentum and conservation of 
energy. In the gaseous phase the state variables P, p,R and T are related 
by the equation of state while for the particulate phase the equations are 
for the particle drag, particle heat balance and the particle equation of 
state. Development of these equations is based on the following assump- 
tions; 

1. The particles are spherical in shape. 

2. The particle internal temperature is uniform. 

3. The gas and particles exchange thermal energy by convection 
and radiation (optional). 

4. The gas obeys the perfect gas law and is either frozen and/ or 
in chemical equilibrium, or is in chemical non-equilibri\am. 

5. The pressure of the gas and the drag of the particles contribute 
to the force acting on the control volume. 

6. The gas is inviscid except for the drag it exerts on the particles. 

7. There are no particle interactions. 
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8. The volume occupied by the particles is negligible. 

9- There is no mass exchange between the phases, 

10. A discrete number of particles, each of different size or chem- 
ical species is chosen to represent the actual continuous particle 
distribution. 

11. The particles ai*e inert. 

The supersonic two-phase solution accepts the starting line provided 
by the internally calculated transonic solution as well as other pertinent 
data supplied through the read function. The equations of motion under the 
assumptions just listed are hyperbolic and permit the use of a forward 
marching scheme; a streamline/normal grid structure is employed where 
the step lengths in the axial and radial directions are under program con- 
trol. Both BCD (printer) and unformatted binary output tapes are pro- 
duced. A Prandtl-Meyer expansion of the gas phase and a free boundary 
calculation are employed to treat the plume flow solution. The run is 
terminated when prespecified problem limits are reached. 

The two -phase flow analysis will treat an extremely wide range of 
operating conditions. With few exceptions the limitations are imposed by 
the theory rather than numerical considerations. In this discussion 
dimension statement sizes which are arbitrarily set are not considered 
a limitation. The true limitations are: 

o Supersonic regions influenced by embedded subsonic regions. 

o Vacuum or limiting expansion limitation — a small region of 
the expansion fan for a vacuum expansion cannot be treated 
where the Mach number is so large that treatment by con- 
tinuous flow assumptions becomes meaningless (this limita- 
tion is both numerical and theoretical). 

o For two-phase flow the lower boundary can only be horizontal 
(i.e., nozzle centerline). 

A complete derivation of the governing equations are available in 
Refs. 1 and Z. 
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A free molecular flow calculation has been provided as an option 
which permits treatment of tbe rarefied regions of the plume. As the gas 
expands it first freezes out the vibrational and rotational modes. During 
this transition the characteristic equations continue to be employed but the 
equation of state is modified. At translational freezing, however, the solu- 
tion switches to an effective source solution. The stream Lines are con- 
sidered straight and the velocity constant. Conservation of mass then 
determines the density while other properties are found from the equation 
of state. This option should not be exercised at this time as it has not 
been checked out with the new modifications which have been added to 
the code. 

2.1 MODIFICATIONS TO THE RAMP CODE 

This section deals with the modifications which were incorporated into 
the RAMP code under this contract. 

2.1.1 JANNAF Standard Plume Flowfield (SPF) Code Interface 

To perform a plume calculation, the SPF code (Ref. 3) requires nozzle 
exit properties as initial conditions. The RAMP code was modified to punch 
or put on tape exit plane data in the format that the SPF code uses. The 
code will punch data for both single and two-phase cases. 

At present, the SPF code can accept a single chemical species distri- 
bution which must be applied at all points across the exit. The RAMP code 
can have a distribution of species, so it was necessary to mass flow average 
the species distributions across the exit plane to produce a single set of 
species. This restriction can be easily removed whenever SPF has the 
capability to handle spatial variations in chemical species. 

The code is set up to punch the species information for the various 
chemical systems that are available in the SPF code. The user must 
specify which system to punch. If the thermodynamic data is input to the 
RAMP code using cards then the species mole fractions that are to be 
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punched for the SPF code must be input for the appropriate system in the 
same order as listed in Table 1 . For finite rate cases these species are 
determined inside the code and need not be input. The user may also input 
another system. Additionally, ideal gas cases can also be handled. 

Table 1 

SPF CHEMICAL SYSTEMS 


System 1 — Hydrogen/Oxygen 
H, H^O, N^, O, OH, O2 
System 2 — Carbon/ Hydrogen/Oxygen 

CO, CO^, H, H2, H2O, N2, O, OH, O2 
System 3 — Carbon/ Hydrogen/Oxygen/Chlorine 

Af 203 (S), CO, CO2, Ci, 0 ^ 2 . H, H2, H2O, HCf, N2, O, OH, O2 
System 4 — Carbon/Hydrogen/ Oxygen/ Chlorine/Fluorine 

CO, CO2, Cf, 0 ^ 2 , F, F2, H, H2, H2O, HCf, HF, N2, 0 , OH, ©2 
System 5 — Hydrogen/Oxygen/Boron 

BO, BO2, B2O3, H, H2, H2O, HBO2, N2, O, OH, O2 
System 6 ~ Hydrogen/Oxygen/Boron/Chlorine/Fluorine 

BF, BF2, BF3, BO, BOCf, BOF, BO2, B2O3, CO, CO2, Ci, CJl^ 

F, F2, H, H2, H2O, HBO2, HCi, HF, 1 ^ 2 , 0 , OH, O2 

2 . 1.2 Shock Capturing Finite Difference Operator 

The original RAMP code has the logic for computing any number of 
right- or left-running shocks using shock fitting techniques. This logic was 
partially checked out under previous efforts. It is desirable to have the 
capability to treat shocks that can occur in some nozzles. Shock capturing 
schemes require no special equations or program logic and are reliable. 
For most nozzle flows, existing shock capturing techniques are sufficiently 
powerful to treat the shocks. Additionally, in order to directly interface 
with the SPF code, it would be desirable to use a shock capturing numerical 
operator. 
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For the above reasons a shock capturing algorithm was added to the 
code. The methodology, equations and grid system which was incorporated 
into tbe code is identical to the SPF code {Ref. 3). To obtain more informa- 
tion on tlie scheme tlie reader is referred to Ref. 3. 

For the shock capturing option, the computational grid are vertical 
lines and a radial coordinate system normalized by the local wall coordi- 
nate. This approach was taken to alleviate the problem of inserting a 
point per line for a Cartesian system. The computational procedure is as 
follows: the forward marching algorithm is used to compute the flowfield 
properties at all interior points. Flowfield properties at the axis and wall 
are obtained via the method-of-characteristics solution. Once the wall and 
axis points have been computed, the coordinate system is transformed and 
the interior points are solved. 

The shock capturing routines have been checked out for ideal gas 
single-phase solutions. Real gas or two-phase flows have not been checked 
out. It is therefore recommended that shock capturing not be used until 
the author has fully verified the logic. When this option has been verified 
the changes to the code will be forwarded to the users. For the present 
time, the regular streamline -normal computation scheme should be used. 

2.1,3 Two-phase, Equilibrium/Frozen, Boundary Layer Analysis 

Nozzle boundary layers are known to influence certain regions of 
nozzle flow and high altitude exhaust plumes (Refs. 4,5 and 6). For nozzle 
designers, the nozzle boundary layer is important in determining the 
thermal loads to the nozzle, performance losses due to heat transfer and 
effects on the nozzle pressure distribution due to the displacement thick- 
ness effect on the inviscid flow structure. For spacecraft designers the 
nozzle boundary layer is important because of its effect on the exhaust 
plume. At high altitudes the nozzle boundary layer causes the plume to 
expand to large angles (approximately 180 deg). In these backflow regions, 
spacecraft and sensitive surfaces are subjected to unwanted contamination. 
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forces, moments and heating rates. For some applications the radiative 
properties of the expanded boundary layer flow is important. Thus, it is 
easy to see that for many applications the nozzle wall boundary layer is an 
important factor. 

To date there have been numerous techniques used to account for 
boundary layer effects of the nozzle and exhaust plume flow fields. Lockheed- 
Huntsviile has an option in the MOC and PIP code (Ref. 7, a version of RAMP) 
which superimposes a power law boundary layer profile based on a flat plate 
boundary layer thickness on the exit plane start line. The profile can contain 
a total temperature variation which has been found to be important (Ref. 7) in 
the back flow region of the plume. This is a very simplified approach and 
the method of determining the boundary layer thickness is very crude. 

Lockheed -Huntsville personnel and other investigators have used various 
boundary layer codes and the inviscid nozzle results to generate a better 
set of boundary layer properties which were then superimposed on the nozzle 
exit plane inviscid results. Next, some investigators have gone a step further 
and iteratively solved the inviscid nozzle and boundary layer to determine the 
effect of the boundary layer on the nozzle flow field. Finally, it is possible 
to do a fully viscous nozzle solution with a method such as Lockheed-Huntsville' s 
GIM code (Ref. 8). However, a fully viscous nozzle solution is beyond the scope 
of this effort. It should be noted that all techniques except the first require a 
tremendous amount of effort and hand calculations to finally obtain a plume 
solution. The purpose of this task was to automate the boundary layer/nozzle 
solution to arrive at an exit plane start line which includes all the effects 
mentioned in tire above. 

The BLIMPJ boundary layer code (Ref. 9) was chosen for the solution 
of the nozzle wall boundary layer. BLIMPJ is a JANNAF standard boundary 
layer code for determining boundary layer effects on the performance of a 
rocket engine. The BLIMPJ code is set up to handle nozzle flows, equilibrium 
or frozen chemistry, uses the JANNAF standard thermochemistry curve fits 
and has numerous ways of treating the nozzle wall that allow it to treat liquid 
engines as well as solid motors with and without ablative walls. At present, 
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the interface between RAMP and BLIMPJ allows only the following wall 
options: adiabatic wall, steady state energy balance or a specified tem- 
perature distribution. The adiabatic wall option is the default option in 
the code. The remaining wall options in the BLIMPJ code could be used 
with a few modifications to the RAMP code. In addition to wall boundary 
condition options, the user may also select laminar, turbulent, frozen, or 
equilibrium flows or let the code default to turbulent, equilibrium. A ' 
turbulent boundary layer is used when the Reynolds number based on mo- 
mentum thickness exceeds 250. 

All input data for the boundary layer solution are generated inside the 
RAMP code except the nozzle wail temperature distribution, if that option 
is specified, Chemical species must be input if the thermochemistry data 
for the inviscid solution was input via cards. 

The BLrIMPJ code was too large for conversion to a module for the 
RAMP code, so if it is desired to generate a nozzle boundary layer and 
obtain a viscous exit plane start line it is first necessary to run the RAMP 
code to calculate the inviscid nozzle solution. The RAMP code generates 
a file with the BLIMPJ input data on it. The BLIMPJ code is then executed 
(no data cards). The BLIMPJ results are printed out and stored on a file 
for subsequent use by the RAMP code. The RAMP code is then reexecuted 
with but a single input card and an exit plane vertical or normal start line 
is generated. The exit plane start line has the results of the inviscid 
nozzle and boundary layer solutions merged together. This capability is 
a considerable improvement over previous methods which required many 
hand calculations, cross plotting, etc. 

For most rocket motors the boundary layer is fairly thin (approxi- 
mately 5 to 10% of exit diameter) and the resultant effect on inviscid flow 
properties is minimal. For these cases one pass through the inviscid 
nozzle solution and boundary layer calculation is adequate. The boundary 
layer results will then be superimposed on the inviscid nozzle solution and 
an exit plane start line with boundary layer effects will he generated which 
can he used to perform a plume expansion. 
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Some low pressure or low thrust motors have boundary layers which 
contain a significant portion of the total mass flow of the system. In these 
cases the entire solution should be iterated for by making two passes through 
the inviscid nozzle and boundary layer calculations. After the initial nozzle 
and boundary layer solution has been completed the nozzle solutions will be 
calculated with the actual wall contour adjusted by the local displacement 
thickness which was determined from the boundary layer calculation. The 
boundary layer solution will then be rerun using the new edge conditions 
from the second nozzle calculation. The results will then be superimposed 
on the original nozzle contour and second flowfield solution and an exit start - 
line will be output or saved. This capability is not presently operational. 
However, for these cases the boundary layer solution is solved a second 
time with new edge conditions taken from the inviscid results at the edge 
of the first BLiIMPJ boundary layer. Thus, the inviscid and viscous results 
match very well at the boundary layer edge. This option is controlled in- 
ternal to the program and requires no user interface, 

Z.1.4 particle Tracing Through Boundary Layer 

Any particulate matter unburned propellant droplets 

in liquid motors) which might enter the nozzle wall boundary layer could be 
influenced by the boundary layer to be expanded into the backflow region. 
Thus, there may be certain motor/nozzle configurations in which it is 
necessary to track particles through the boundary layer in order to obtain 
the most*accurate representation of the nozzle flow properties at the exit 
plane. 


Lockheed-Huntsville added a particle streamline tracing module into 
the nozzle code. The code which was used as the basic building block of 
this module was used in predicting the lUS, SSUS particle distributions 
published in Ref, 10. This code uses initial particle properties (velocity, 
flow angle temperature and density) and traces the particles through a 
known flow field. This option is user-selectable. The basic method 
used is shown in Fig, 1. 


10 


LOCKHEED-HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D784753 



i 


Fig. ^ “ Schematic of Particle Trajectory Tracing Method to be Used 
in the Inviscid Nozzle Code 
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After the inviscid nozzle solution and boundary layer solutions have 
been completed, the RAMP code is re-executed and logic is handed off to 
the particle trajectory tracing module. The first step is to search the 
boundary layer edge location in the inviscid nozzle solution to determine if 
any particles are found to penetrate the boundary layer. Then for each par- 
ticular size group which has been found to penetrate the boundary layer a 
distribution of individual particle properties will be determined (interpolated 
off nozzle solution tape) at several stations along the boundary layer edge. 
This will be done for each size group. Next a map of the boundary layer 
properties will be generated using the boundary layer data tape. Finally, 
each particle streamline will be traced through this mapped flow field and 
particle properties will be determined at the previously determined exit 
plane. The effect of the particles on the boundary layer flow properties 
due to energy and momentum interchange are not considered. 

Finally, the boundary layer results, the inviscid nozzle solution and 
the particle properties in the boundary layer are merged, and an exit plane 
start line for a RAMP restart or SPF run will be generated. 

2.1.5 Boundary Layer Expansion at Lip 

Another required capability of the nozzle code is a fully automated 
technique to expand the boundary layer at the lip. After an exit plane start 
line has been generated it is very straightforward to use the Prandtl-Meyer 
relationship to expand the flow. The key point here is to get a fully super- 
sonic start line. Numerous investigators have made various assumptions 
to obtain start lines for plume solutions that account for the nozzle wall 
boundary layer. The various techniques will be discussed briefly. 

Perhaps the most simple method of obtaining a supersonic start line 
is to discard the subsonic portion of the boundary layer (Refs. 4 and 11) and 
adjust the spatial distribution of boundary layer points up to the lip. This 
method throws mass away from the system and will affect the results in the 
backflow region. For nozzles which have thin turbulent boundary layers, 
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these effects are probably minimal. For very viscous nozzles (low pres- 
sure) or laminar boundary layers the effects of this assumption could 
severely compromise the plume results in the backflow region (greater 
than 90 deg). It should be noted that this technique could consider total 
temperature variations through the boundary layer, 

Seubold (Ref. 12) used the technique of replacing the subsonic portion 
of the boundary layer with mass average properties at a slightly supersonic 
value. This resulted in a layer of constant properties near the wall. The 
resultant plume flow had the correct mass flow but did not consider the total 
temperature variations. 

The first two techniques assumed that the static pressure through the 
boundary layer is constant. Near the lip however the subsonic portion of 
the boundary layer is influenced by the expansion process (Refs. 13, 14 and 
15) as shown in Fig. 2. For highly underexpanded flows the sonic line has 
been, found to attach to the lip (Ref. 13) so that the flow at the wall must 
rapidly accelerate when it gets near the lip and the static pressure rapidly 
decreases. For overexpanded flows the subsonic flow merely stays subsonic 
as it negotiates the lip so that downstream flow conditions can feed back up 
into the boundary layer. 



Fig. 2 - Mach Number Contours for an Underexpanded Axially 
Symmetric Nozzle (Case 6) 
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Cooper (Ref. 16) went a little further than the first two approaches as 
far as treatment of the subsonic flow expanding to the sonic lip conditions. 

He assumed that the flow (M=0 at wall) just upstream of the lip expanded to 
M = 1 at the lip. Since the static pressure upstream of the lip at the wall 
was the same as the wall total pressure and the freestream static pressure, 
the only way for M = 1 to be reached at the lip was for the wall static pres- 
sure to drop. Therefore, there will be a static pressure distribution across 
the flow at the lip. It was then assumed that the static pressure at specified 
point in the boundary layer is not affected by the expansion and the pressure 
distribution is faired into the known lip conditions. This method is not totally 
rigoi'ous but does include boundary layer effects and has been, shown to result 
in good comparisons with some experimental data in the backflow region. 

Finally, there are exact solutions to the corner flow problem. Bird 
(Ref. 13) uses a direct simulation Monte Carlo method. He set his model up 
so that it would compute the entire nozzle or the region near the throat. He 
predicts the attachment of the sonic line for underexpanded nozzles. Baum 
(Ref. 15) uses a finite difference method along with boundary layer equations 
to describe the subsonic portion of the boundary layer. His application was 
for base flow about blunt bodies and compares well with data. Finally, the 
previously mentioned GIM code can be used to exactly solve the corner prob- 
lem. All three of these methods are very complex and are outside the scope 
of this effort. 

In the present version of the code the method of Seubold was incorporated. 
Total temperature variations are included and the mass flow in the boundary 
layer is conserved. Plume calculations using this method were made for a 
5 Ibf bipropellant motor. The mass flux predictions compare very well with 
experimental data. 

It is anticipated that further effort will be spent on the lip region so 
that better calculations may be obtained in the backflow region of the plume. 
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2.1.6 Variable Oxidizer -to-Fuel Ratio Transonic Solution 

Solution of the subsonic transonic region of a liquid rocket engine can 
vary in complexity from a simple one -dimensional variable O/F streamtube 
analysis (Ref. 17) to the most detailed model such as that of the Distributed 
Energy Release (DER, Ref. 18) model. The streamtube analysis performs a 
multizone one -dimensional calculation to the sonic point given a known O/F 
distribution just downstream of the injection face. The DER program is a 
complex model which is Initiated upstream of the injector face and continues 
the solution up through the sonic line. The DER code was used to initiate 
nozzle solutions in Ref. 19 but is not particularly easy to use and input and 
requires a good bit of experience of the user to successfully execute. For 
these reasons it will not be utilized in the nozzle code. On the other hand, 
a one “dimensional streamtube analysis does not account for two-dimensional 
effects. A time -dependent scheme (Ref. 20) is a compromise between these 
two schemes. The approach includes the radial momentum equation which 
results in a set of mixed partial differential equations. The solution pro- 
cedure is an unsteady time-dependent finite difference technique with equi- 
librium chemistry. This technique has both the equilibrium and variable 
O/F chemistry option and has been utilized (Refs. 21 zuid 22) previously with 
excellent results. This code has been incorporated into the RAMP code as 
a module and has been executed for all combinations of ideal and equilibrium 
chemistry, constant and variable O/F distributions. Additional inputs re- 
quired by this module are; (1) combustion chamber geometry down to the 
throat; (2) average O/F ratio, (3) location of the initial data surface to start 
the transonic solution; (4) location of the nozzle throat, and; (5) the O/F dis- 
tribution at the initial data surface. 

2.1.7 Improved Two-Phase Transonic Solution 

The original transonic module which was incorporated into the RAMP 
code could handle throat radii of curvature to throat radius ratios above 1.5. 
Many solid motors have radii of curvature ratios smaller than 1. To alleviate 
this limitation the improved approximate transonic module was taken from the 
new Standard performance Prediction program (SPP, Ref. 23) and put into the 
RAMP code. 
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Many transonic solutions using this module have been run for varying 
upstream and downstream radii of curvature ratios. Unfortunately a prob- 
lem was encountered for the smaller radii of curvature cases. It was not 
possible to obtain a true normal for starting the supersonic solution within 
the domain of the transonic solution. For this reason the old module was 
put back in the code and the limitation of radii of curvature ratios greater 
than 1.5 must remain in effect. Work is continuing to interface the transonic 
with the supersonic solution and the new modules will be forwarded to the 
users as soon as they are checked out. 
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2. 2 CAPABILITIES AND LIMITATIONS 

The RAMP computer program described in this document can be used 
to solve a wide variety of problems associated with real gas, supersonic, 
compressible flow. Some of the more important, basic capabilities of the 
existing program are outlined below: 


• The gas may be ideal or real. If the gas is real; frozen, 
equilibrium, or non- equilibrium chemistry assumptions 
can be made. The effects of oxidizer/fuel gradients may 
be considered. 

• Two-dimensional or axisymmetric flow problem geometries 
can be used. 

• Both upper and lower boundaries can be solid or free. {A 
solid boundary can be approximated by either a conic or 
polynomial equation.) (Two -phase problems require the 
nozzle centerline as a lower boundary). 

• A nozzle wall may be curve fit with discrete points. 

• Any number of expansion corners can be considered on 
either the upper or lower wall. 

• Various methods for obtaining an initial start line are 
utilized. 

1. The program will calculate a one -dimensional 
start line anywhere in the nozzle. 

2. The program will calculate a start line at points 
within the nozzle necessary to conserve mass. 

3. Data on a normal surface can be input at points 
across the flow field within the nozzle or in the 
plume. 

4. An exit plane startline can be punched. 

5. The program can be restarted from the startline 
punched in 4 above (except viscous exit plane) 

6. Single phase start line can be calculated using 
a transonic solution which can handle ideal gas, 
equilibrium chemistry for both constant and vari- 
able oxidizer/fuel ratios 

7. Two phase start line can be set up using a tran- 
sonic solution. 
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• Hypersonic or quiescent approach flow options may be 
used. 

• Exit to ambient pressure ratios from over-expanded to 
highly under -expanded are possible. 

• Displacement of the axis of symmetry from the center 
of flow {i, e., the plug nozzle flow field) is possible for 
gas only cases). 

• There is presently a maximum of 100 points on a normal 
and 50 input points. 

9 Reacting gas solutions which are in chemical equilibrium 
have been facilitated by modifying the TRAN72 (Refs. 24, 

25) computer program as described in Section 2 of Ref. 2 
to provide binary tape and punched output of its equilibrium 
or frozen real gas calculations at any desired O/F ratio(s) 
or total enthalpy(s). The RAMP program has the capability 
for selecting the proper case from a large set of real gas 
properties cases stored on a master tape. The method of 
generating this master tape is outlined in Table 2. Cases 
stored are uniquely identified by some characteristic of 
the particular gas under consideration. For example, a 
LOX/IjH^ system may be identified by the following: 

Gas Type Mixture Ratio Chamber Pressure 
02/H2 O/F = 1.5 - 8.0 PC = 546.0 

New cases of general interest may be added to the.master 
tape; however, ad hoc cases should be prepared on a sepa- 
rate tape. Tape preparation sequence and communication 
with the RAMP program is diagrammed in Table 2. 

e Once the gas -particle flowfield solution has been obtained, 
the output tape may be used by the RAMP Radial Hookup 
Program (described in Appendix A of Ref. 2) which deter- 
mines the radial variations of flowfield properties across 
the nozzle and plume flow fields at constant axial stations. 
The Plume Impingement Program (PLIMP) (Ref. 26) may 
also be run to determine the effects of the rocket exhaust 
plume on objects immersed in the plume. Sequencing and 
communication of auxiliary programs with the RAMP pro- 
gram is shown in Table 2. 

o Two-dimensional or axLsymmetric solutions are selected 
by simply loading a control word in the program input 
data. This integer (0 or 1) is then multiplied by the term 
containing (l/r) in the governing differential equation. By 
appropriate description of the flow boundaries, it is possible 
to change from a solid to free boundary on either the upper 
or lower walls. Conversely, it is not possible to change 
from a free to a solid boundary on either wail. 
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, Table 2 

SEQUENCING AND COMMUNICATION OF AUXILIARY PROGRAMS WITH THE RAMP PROGRAM 
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• A real gas nozzle boundary layer solution can be 
performed with no interface between the user and 
the RAMP code results. A frozen, equilibrium 
chemistry, turbulent or Laminar solution is possible. 

• The effect of a nozzle wall boundary layer on par- 
ticles which enter the boundary layer is treated. 

® An exit plane start line can be generated for the SPF 
code which includes two phase flow and boundary 
layer effects. Additionally exit plane data is avail- 
able for use in other codes. 

« For two phase transonic solutions the nozzle throat 
radius of curvature/throat radius ratio must be the 
same upstream and downstream of the throat and 
greater than 1.5. 

• A restart at the exit plane using a start line that 
includes the nozzle boundary layer merged with 
the inviscid results is not possible. Additional 
logic for mesh and step size control needs to be 
added. 

• Multiple passes through the inviscid nozzle solution 
with the nozzle contour modified by the BLiIMPJ solu- 
tion displacement thickness is not possible. 

« Shock capturing numerical analog logic is not fully 
checked out and should not be used at this time. 


2.3 USER'S INPUT GUIDE FOR THE RAMP PROGRAM 

This section outlines in detail the procedures for using the Reacting 
and Multi-phase (RAMP) Computer Program. Each card and its use is 
explained in Section 2.3,1. The program magnetic tape assignments are 
given in Table 3. 

2.3.1 RAMP Program Input Information 

The input data are organized into sections determined by their use. 
The description of these cards is given below. 
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RAMP Computer Program Input Instructions 


Card 1 

Column 

5 


Overall Run Control Card 
Parameter Value 

ITRCE 0 


1 


10 ISTART 0 


1 

Z 

3 


Format (l6l5) 

Description 

All input data are coming 
from cards except as noted 
for rCON(l) and ICON(2). 

Input cards will be written 
on unit 13 for future use 

RAMP code has already 
been executed for a nozzle 
and all input data from noz- 
zle run is on unit 13. Nozzle 
boundary layer has been com- 
pleted, and the code will gen- 
erate an exit plane start line 
which includes boundary layer 
and particles in the boundary 
layer (if present). 

Do not generate an exit plane 
start line that includes the 
boundary layer. Inviscid start 
line is controlled by Card 5. 

Punch a viscous SPF start 
line. 

Punch a viscous normal start 
line at exit. 

Punch a viscous normal and 
SPF start line. 


4 


Cards 2-4 

Column 

1-240 


Problem Description 
Required 

Parameter 

HEADER 


Card 5 


Run Control Card 
Required 


Punch a viscous normal (or 
vertical if ICAPTR>0) and 
SPF (if ISKPY>0) start line 
and restart at exit. Note; 
This option is not operational. 

Format 3(20A4) 


Description 

problem description may be 
put on three cards; however 
only the first 120 columns will 
be printed while all 240 char- 
acters will be written on the 
data tape. All three cards 
must be present even if blank. 

Format l6l5 (Right Adjusted) 
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Card 5 (Cont'd) 

Column Parameter Value 

5 ICON(I) ■ 1 

Gaseous thermo- 
dynamic data con- 
trol parameter 


2 


3 


4 


8-9 NTAPE 


N 


10 ICON{2) 0 

Start line control 
parameter for gas 
only solution 

2 

or 3 


4 


Startline control 0 

parameter for gas- 
particle solution ^ 


Description 

The gas composition is either 
chemically frozen and/or in 
chemical equilibrium. The gas 
properties are read directly 
from cards 9, 10, 11 and 12. 

Same as ICON(l)=l except gas 
properties are read directly 
from a data tape mounted on 
FORTRAN unit 10. 

The gas composition is in chem- 
ical non- equilibrium. The gas 
properties are determined, as 
a function of temperature in 
thermodynamic data tables in- 
put on card 14. 

Same as ICON(l)=3 except gas 
composition is chemically 
frozen. 

If ICON(2)=2, tape unit numiber 
for startline if not input from 
cards. The program defaults 
to unit 5 (read cards) for 
ICON{2)=2. If ICON{2)/2 and 
a two-phase transonic solution 
is being performed NTAPE is 
the unit on which the ti'ansonic 
startline will be written. In this 
case the program defaults to 
unit 8 

Generate straight startline with 
Mach number given. 

Generate source startline with 
a/ a given 

Startline input from cards or 
tape. 

Generate startline by conserva- 
tion of mass using a linear Mach 
number distribution. 

Generate gaseous startline using 
variable O/F transonic module. 
This option requires cards 19c, d,e 

Generate startline using tran- 
sonic approximation. 

Startline input from cards. 
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C ard 5 (Conf d) 

'l3 ICON(3) 

Control manner in 
which points along 
the startline are 
s paced 


14-15 ICON(3) 

16-ZO ICON(4) 

Number of upper 
boundary equations 


0 Points are spaced according to 
a sine distribution. 

1 Points are evenly spaced (recom- 
mended) 

2 Points are evenly spaced on a 
circular arc based on the input 
value of the upper limit of the 
startline (card 19 CORLIP(2)) 

NOTE: This option is necessary 
only if program is to set up its 
own gaseous startline. 

Number of startline points. 
Maximum of 50 (right adjust)”' 

Upper boundary specification 
indicator. If specifying upper 
boundary by equations, set equal 
to number of equations to be used 
Maximum of 100. Right adjust. 


Option for ICON{4) when upper boundary is described by individual points and 
slopes 

ICON(4) INOOO + Number of discrete points (no boundary equation 

following last point) (slope at each point in radians) 

2N000 + Number of discrete points + 1 (an upper boundary 
equation follows last point) (slope at each point in 
radians ^ 


3NOOO + Number of discrete points (no boundary equation 

following last point) (slope at each point in degrees) 

4N000 + Number of discrete points + 1 (an upper boundary 
equation follows last point) (slope at each point in 
degrees) 

N ~ number of points to use for Lagrangian Integration (5 max). 

If N IS set to zero, a linear assumption will be made. 

NOTE: If a nozzle is being run the throat must also be specified by discrete 

points. 


NOTE: If particles are present and supersonic startline is generated by 
transonic approximation then total number of points on startline may be 
adjusted by transonic program depending on particle distributions. 
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Card 5 (Cont'd) 




Column 

Parameter 

Value 

Description 

21-25 

lCON(5) 

Number of lower 
boundary equations 


Lower boundary specification 
indicator. Same description 
and option as ICON(4), 

27-28 

ICON(3) 

N 

N is the boundary equation num- 
ber at which the SKIPPY start 
line is to be punched. The code 
uses the XMAX (see card 7) for 
this boundary equation as the 
station to punch. 

29 

lUNIT 

0 

Punch SPF start line 



1 

Put SPF start line on unit 12. 

30 

ISKPY 

0 

Do not generate SPF start line. 



1 

Chemical system 1 SPF start 
line (see Table 1) 



2 

Chemical system 2 SPF start 
line (see Table 1) 



3 

Chemical system 3 SPF start 
line (see Table 1) 



4 

Chemical system 4 SPF start 
line (see Table 1) 



5 

Chemical system 5 SPF start 
line (see Table 1) 



6 

Chemical system 6 SPF start 
line (see Table 1) 



7 

Input chemical system via card 



8 

Ideal gas SPF start line. 

33 

ICAPTR 

0 

Streamline -normal solution 



1 

Shock capturing numerical 
operator (This option not fully 
Operational. ) 

34 

IBL 

0 

Do not generate data for bound- 
ary layer calculation. 



1 

Generate data for boundary 
layer calculation and store it 
on unit 1. 

35 

ICON(7) 

0 

Two-dimensional flow problem 
geometry. 



1 

Axisymmetric flow problem 
geometry. 
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Card 


{Conf d) 

38 INOZ 2 

39 ICON(8) 0 

Data output con- ^ 

trol, used in con- 
junction with ICON(l6) 

40 ICON (8) I 

2 

3 

42 MORFT 
Compliments I CON (9) q 

1 

2 

3 

4 

43 lCON{9) 0 

Units indicator , 


44-45 ISPECS 
48-50 ICON(IO) 

25 


Calculations terminated at 
nozzle exit. Set if boundary 
layer is to be calculated. 

Full printout 

Print only boundary, shock, input, 
Prandtl-Meyer, and particle limit- 
ing streamline points. 

Print 1 line (R, X, M, 0, S and 
shock angle) 

Print above plus Mach angle, 

P, P, T, V. 

Print all of above plus MWT, 
y, TO*, PO* 

For English system of units. 

Dimensions are in feet. 

Dimensions are in inches. 

For metric system of units. 

Dimensions a‘re in centimeter's. 

Dimensions are in meters 

System of units specified by user. 

Use English system of units. 

Use metric sjrstem of units. 

This option controls the units in 
which the flow field is calculated. 

The program assumes that the 
boundary equations are input in the 
same units as the units indicator 
(ICON(9)). This option will not 
override the units specification on 
cards 9 and 31 but will convert the 
units of the gas and particle thermo- 
dynamics to correspond to the units 
of this indicator. 

Number of discrete particle sizes 
used to represent particle distri- 
bution (10 max). If gaseous only 
flow set equal to 0 (right adjust). 

Maximum iterations allowable 
for each point in flow field. If 
set to 0 program assumes value 
of 100. Right adjust. 
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Column 

Parameter Value 

Description 

51-55 

rcoN(ii) 


Case number printed at top of 
each page. 

60 

ICON(IZ) 

0 

Calculate shock wave. 



1 

No rotation option. 

61 

ICON(13) 

0 

Flowfield data will be output on 
FORTRAN unit 3. Must be set 
if generating boundary layer 



1 

Data will not be written on tape. 

65 

IC0N{13) 

0 

Free molecular calculations will 
not be considered. 



1 

Free molecular calculations 
will be considered. 

68-70 

lCON(l4) 

0 

No intermediate printout in 
solution iteration. 



N 

Pr^t intermediate results for 
line. Right adjust. 

71-75 

lCON(15) 

0 

No intermediate printout. 



M 

Print intermediate results from 
point on each line from the 
Ntti (ICON (14)) line on. Right 
adjust. 

76 

ICON(l6) 

0 

No punched cards output. 



1 

Punch data line at nozzle exit 

77-78 

lCON(l6) 

0 

Print every line. 



N 

Print every Ime (use with 

ICON (8)). Put 0 in column 77 
if N < 10. 

79-80 

lCON(l6)* 


Time (SEC) before end of allotted 
run time when new startline is to 
be punched. Put 0 in column 79 
if time less than 10 seconds. 

Card 6 

Finite Kate Chemistry 
Run Control Card (Re- 
quired if ICON(l)>2) 


Format 815 (Right Adjusted) 

Column 

Parameter Value 

Description 

1-5 

NT 


Number of temperature points in 
thermodynamic data tables. 

, 6-10 

NS 


Number of gaseous species (ex- 
cluding 3rd bodies) 

* 

NOTE: Applicable for Univac 1108 onlyr. 
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Column 

Parameter 

Value 

Description 

11-15 

NM 


Number of 3rd bodies. 

16-20 

NR 


Number of reactions specified. 

25 

NPRINT 

0 

No intermediate printout in 
chemistry calculations. 



1 

Echo print of input data. 



2 

Print intermediate results of 
chemistry calculations. 


30 ICTAPE 0 Species concentrations for start - 

line read directly from cards. 

1 Species concentrations read di- 

rectly from a data tape mounted 
on FORTRAN unit 10. 

35 KGUP >2 Number of normals calculated 

before finite rate chemistry 
contributes to dS and dH. 

40 IDIDO 0 Uniform species concentrations 

along startline. 

1 Non-uniform species concentra- 

tions along startline. 


Card 7 Upper Boundary Description 

Required* 

If ICON (4) < 10000 use following format (II, 3X,I1, 5X, 6E10.6). 

Column Parameter Value Description 

1 IWALL{K, 2) 1 Conic equation 

R = A[(B+CX+DX^) ' +e] 

Represents throat region. (See 
page 3-27 for an example and 
description.) 

2 Polynomial equation 

R = AX'^+BX^+CX^+DX+E 

3 Free boundary equation 

P = P (1+E X)(l +7 (M sin(6^-0 ))^) 
00 00 'oo 00 B 00 

(See page 3-27 for an example and 
description.) 


6 Same as IW'ALL=3 except oblique 

shock solution for plume boundary. 
Use if 1.5 < M < 5.5. 

00 

^ transonic calculation is to be performed the boundary 

conditions upstream of the throat have to be input. 


27 



LMSC-HREC TR D7 84753 


Column Parameter Value 

5 ITRAN(K, 2) 0 

1 

2 

11-20 WAL,LCO(K, 1, 2) 

21-30 WALL,CO(K, 2, 2) 

31-40 WALLCO(K, 3, 2) 

41-50 WAL,LCO(K,4, 2) 

51-60 WALLCO(K, 5, 2) 

61-70 WALLCO(K, 6, 2) 

71-80 RSTAR 


Description 

No discontinuity follows this 
equation . 

Expansion corner follows. 
Compression corner follows. 

Coefficient A or P {psfa or 
2 00 
N/m ), (Units must be con- 
sistent with R in ft or m.) 

B or y 
oo 

C or M 

00 

D or 9 (deg) 

E or E 

00 

Maximum value of X applicable 
to equation (feet if ICON (9) =0 
meters if lCON(9)=l). 

Throat radius (ft or m) required 
only on card for last equation. Re- 
quired for 2 -phase trans. solution 
and boundary layer solution. 


If 10000 < ICON (4) < 20000 use following format (15, 5X. 3E10.6,I5, 
5X, 3E10.6). 


Column 

Parameter 

Value Description 

5 

1 TRANS (K, 2) 

Same as before. 

11-20 

WALLCO(K, 3, 2) 

Axial displacement (X) of point 
K (ft or m). 

21-30 

WALLCO(K, 1, 2) 

Radial displacement (R) of point 
K (ft or m) , 

31-40 

WALLCO(K, 2, 2) 

Wall angle (9) at point K (rad). 

45 

ITRANS(K+1, 2) 

Same as before. 

51-60 

WALLCO(K + l, 3, 2) 

X at point K-rl (ft or m) . 

61-70 

WALLCO{K + l, 1, 2) 

R at point K+1 (ft or m). 

71-80 

WALLCO(K+l, 2, 2) 

6 at point K+1 (rad). 


NOTE: Card 6, in the above format, is repeated for each equation until all 
necessary equations have been input. That is, repeat Card 7, in succession 
in order of increasing XMAX, for K=l, 2, . . . ,ICON(4). All units for lengths 
for two-phase calculations are consistent with ICON{9), otherwise units for 
lengths are input at user's discretion. 
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Repeat Card 7j in above format, m succession, and in order of increasing 
X, until all required points have been input. 

If 20000 < ICON(4) < 30000 the above format is used except the last 
segment of the upper boundary is input via an equation. The equation is input 
with the format for ICON (4) < 10000 except the throat radius RSTAR is not 
required. 

If 40000 < ICON{4) < 50000, the format for 20000 < ICON(4) < 30000 is 
used except 9 is input in dimension of degrees. 

) 

Card 7a Format E10.6 

Card 7a IS used only when running a two -phase case where the upper 
boundary (nozzle wall) is specified by discrete points (ICON (4) > 10000). Do 
not input this card for any other cases. 

Column Parameter Description 

1-10 RSTAR Throat radius (ft or m) . 

Card 8 Lower Boundary Description 

Required 

The formats and options for Card 8 are controlled by ICON (5) and are 
the same as for Card 7 (Upper Boundary) with the following exceptions; (1) 
the distance from the nozzle throat to the center (XsO) of the coordinate sys- 
tem for the wall equations is read m place of RSTAR. This is only necessary 
for two-phase cases where a transonic solution is desired and where X/^0 at 
the nozzle throat. This distance is positive if the center of the coordinate 
system is downstream of the throat and negative if the center of the coordinate 
system is upstream of the throat. It is not possible to run a two-phase case 
with the lower boundary specified by points, therefore there is no Card 7a; 

(2) the indices of the parameters are (““,~, 1) instead of (~, ~, 2 ), e.g., 
WALLC0(K, 1, 1) instead of WALLCO(K, 1, 2). A nozzle throat region showing 
the coefficients of a circular throat and free boundary are shown in the sketch 
on the following page. 
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RC “ radius of curvature of the 
circular arc of the throat 
RT = throat radius 
XO ^ axial distance from the origin 
of the coordinate system to the 
throat 

0 == throat divergence angle corre- 
sponding to the maximum value 
for which the throat conic equa- 
tion applies 


The conic equation for this case would have the following form: 

A = -1 for an upper equation, +1 for a lower equation (~1 for this case 
B = RC^ - XO^ 

C = 2X0 
D = -1 

E = -(RC + RT) 

Xmax = RC sinQ + XO 


An example of a free boundary is shown in the sketch below. 



The free stream approach flow is 
inclined at 15 deg to the plume with 
a gamma (7) of 1.4, a Mach number 
of 10, and a static pressure of 0.1 
psta. 


PINF = 0.1 (psfa) 

E = 0 (No pressure variation with axial distance) 
GAMMAINF = 1.4 
MINF = 10 
THETAINF = -15° 
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Card 9 Gas Property Control Format 6A4, 5X, A3, 6X, 12, 3X, IZ 

This card is required whether gas data input by cards or tape. 


Column Parameter 

1-24 ALPHA (I) 

30-32 UNITS 

{Independent 
of ICON (9)) 

39-40 lOF 


44-45 IS 


Description 

Gas name for real gas on tape. If 
inputting gas data via cards, may 
be any name. 

ENG Input gas data with English 
units (cards only). 

MKS Metric units {cards or tape). ■ 

Number of O/F tables for gaseous 
only solution or 

number of gas total enthalpy tables , 
for two -phase solution 

Number of entropy tables per lOF 
entry, 1 for gas, 2 maximum for gas 
chemical equilibrium solution. 


Card 10 Mixture Ratio or Total Format E10.6, 8X, 12 

Enthalpy {This card is 
not used if ICON{l)>2) 


Column Parameter 

1-10 OF RAT (M) 


11-20 K0(M) 


Card 11 Entropy 

{This card is not 
used if ICON{l)>2) 

Column Parameter 

1-10 STAB {M, I) 

19-20 IVTAB{M, I) 


Description 

For gaseous only flow input O/F ratio, 
for particle flow input gas total enthalpy 
(cal/gm for metric, Btu/lbm for English, 
units specified by Card 9). 

For gas only transonic case input 
total enthalpy {cai/gm for metric, 
Btu/lbm for English). If real gas 
get data from NASA-Lewis run. 

If ideal HO = C T where C = 

YR/iY-1) ^ ° P 

Format ElO.6, 8X,I2 


Description 

Entr,opy of gas {cal/gm-®K or Btu/lbm- 
°R, units specified by Card 9). 

Number of Mach numbers for this 
entropy value {13 max). 
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12 

Gas Properties 
(This card is not used if 
ICON(l)>2; units specified 
by Card 9) 

Format 8E10.6 

Column 

Parameter 

Description 

1-10 

XSIDIM(l) 

Mach number associated with 
above entropy. 

11-20 

XSIDIM(2) 

Molecular weight of gas (gm/g- 
mole or Ibm/lb-mole). 

21-30 

XSIDIM(3) 

Gamma (C /C ) 

31-40 

XSIDIM(4) 

Temperature (°K or °R) 

41-50 

XSIDIM(5) 

Pressure (atm) 

51-60 

XSIDIM(6) 

Prandtl number (dimensionless) 

61-70 

XSIDIM(7) 

Absolute viscosity (poise) 

71-80 

XSIDIM(8) 

Ideal gas (1 velocity cut per 
table) — viscosity temperature 
exponent. Real gas — C^ (cal/ 

gm-°K or Btu/lbm-°R). 


To illustrate the arrangement of Cards 10, 11 and 12, let IOF=2 and IS=2; 
then the proper arrangement is: 

Card 10 
11 

12 (1-13 such cards) 

11 

12 (1-13 such cards) 

10 

11 

12 (1-13 such cards) 

11 

12 (1-13 such cards) 


Card 13 Gas Properties Format 3E10.6 

(This card is required 
if ICON(l)>2) 

C olumn Parameter Default Description 

Value 

I- 10 PR 0.7 Prandtl number (dimensionless). 

II- 20 VISO l.OE-04 Absolute Viscosity (poise). 
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Card 13 (Continued) 

Column Parameter Deiault Description 

Value 

21-30 EX 0.6 Viscosity temperature exponent. 

Card 14 Gas Thermodynamic Data 

(The following cards are 
required if ICON(l)>2) 

The following set of cards contain species thermodynamic data. The 
first card contains the species name, molecular weight and heat of formation. 
The second and remaming cards contain the temperature and corresponding 
specific heat, entropy and enthalpy for that species. Two temperatures and 
corresponding thermodynamic data are placed on each card. The input table 
can contain up to a maximum of 30 temperature points. The data are input 
exactly as presented in the JANAF tables (Ref. 27) with the temperature points 
being the same for all species. Cards 14.1, 14.2, 14.3, etc., are repeated for 
each species. If a SPF start line is to be punched (0<ISKPY<7) then the species 
must be input in same order as they appear in Table 1. 


Card 

Column 

Description 

F ormat 

14.1 

1-8 

Name of first species (Left adjusted) 

2A4 


11-20 

Molecular weight 

E10.3 - 


21-30 

Heat of formation, h-^. (kcal/mole) 

E10.3 

14.2 

1-10 

First temperature point (°K) 

F10.4 


11-20 

c (cal/mole-°K) 

F10.4 

- 

21-30 

(cal/mole -^K) 

F10.4 


31-40 

h. -h^gg (kcal/mole) 

F10.4 


41-50 

Second temperature point (°K) 

F10.4 


51-60 

c (cal/mole-*^K) 
^ i 

F10.4 


61-70 

1 

(cal/mole-'^K) 

F10.4 


71-80 

h -h__„ (kcal/mole) 

F10.4 

14.3 

1-10 

Third temperature point 

F10.4 


etc. 
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Cards 15 Catalytic Species Weighting Factor Data 

(The following cards are required if 
ICON(l)>2 and NM>0) 

The following set of cards specify the catalytic species (Ml, M2, M3, . . .) 
and their respective composition in terms of the species participating in the 
reactions. Weighting factors must be read in the same order in which the 
thermodynamic data sets are read. 


Card 

Column 

Description 

F ormat 

15.1.1 

1-8 

AID(NS+1) — Name of first catalytic species 
(e.g.. Ml) (Left Adjusted) 

2A4 

15.1.2 

1-5 

WF(1, 1} — Weighting factor of first species 
(for first catalytic species). Set weighting 
factor to zero for any reactant which does 
not contribute to the respective catalytic 
species . 

16F5.2 


6-10 

« 

• 

WF(l, 2) — Weightmg factor of second species 
contributing to first catalytic species. 



75-80 

WF(1, 16) — Weighting factor of l6th species 
contributing to first catalytic species. 


15.1.3 

1-5 

WF(1, 17) — Weighting factor of 17th species 
contributing to first catalytic species, etc. 

16F5.2 

15.2.1 

1-8 

AID(NS+2) — Name of second catalytic 

2A4 

15.2.2 

1-5 

WF (2, 1) — Weightmg factor of first species 
contributing to second catalytic species, etc. 

16F5.2 

16.NM.1 1-8 

AID(NS-fNM) — Name of last catalytic 
species, etc. 

2A4 

Cards 

16 Chemical Reaction Mechanisms 

(The following cards are required 
if ICON(l)>2 and NR>2) 



The following 

set of cards specifies the chemical reaction mechanisms 

for a particular problem, one card for each reaction. No particular order is 

required. Species 

names are Left adjusted. 


Card 

Column 

Description 

F ormat 

16.1 

1-8 

Species A 

2A4 


9 

+ sign 



10-17 

Species B (or M) 

2A4 
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Cards 16 (Continued) 


Card 

Column 

Description 

F ormat 

16.1 

18 

+ sign 



19-22 

Blank (or M) 

A4 


23 

= sign 



24-31 

Species C 

2A4 


32 

+ sign (if needed) 



33-40 

Species D (or M) 

2A4 


41 

+ sign (if needed) 



42-49 

Species E (or M) 

2A4 


50-51 

Reaction type, 1 to 12 

12 


52 

Rate constant type, 1 to 5 

11 


53-59 

A, pre- exponential factor (cm- 
particle-sec units) 

E7.1 


60-64 

N, temperature exponent 

F5.2 


65-74 

B, activation energy (cal/mole) 

FlO.i 


75-80 

M, temperature exponent 

F6.2 

16.2 


Next reaction 


16.NR 


Last reaction 


Cards 

17 Startline Data Format 

7E10.3 


(The following cards are required 
if ICON(l)>2 and ICTAPE=0) 

The following cards contain the species mole fractions on the startline. 
Mole fractions must be read in the same order m which the thermodynamic 
sols are read. 

Card Column Description 

17.1 1-10 Mole fraction of first species at the first point on the 

startline. 

61-70 Mole fraction of seventh species at the first point on 
the startline. 
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Cards 17 (Continued) 

Card Column Description 

17.2 1-10 Mole fraction of eighth species at the first point on the 

startline. 


61-70 Mole fraction of the fourteenth species at the first point 
on the startline. 


etc. 


Cards 17,1 and 17.2, etc., are repeated for each point on the startline. 
l*or a uiuli>rm sfurllino (IDIDO^O), aioLu iraclioiib arc read ior I point only. 


Card 18 Chamber Condition Data Format 2E10.3 

(This card is used if ICON(l)>2 
and ICTAPE=0) 


Column 

1-10 

11-20 


Parameter Description 

PC Chamber pressure (atm) 

TC Chamber temperature (°K) 


Card 19 Startline Data Format 8E10.6 

(This card is not used if ICON (2) =2 
or for gas particle flow). 


Use Card I9a if ICON(l)<2 and ICON(2)<3. Use Card 19b if ICON(l)>2. 
Use Caidb 19c, d, c xi ICON (2)^ 4 


Card 

19a 


Column Parameter 

I- 10 CORLIP(2) 

II- 20 CORLIP(6) 

21-30 CORLIP(4) 

31-40 CORLIP(5) 

41-50 CORLIP(8) 


Description 

Axial coordinate of upper limit of startline 
(ft or m, see Fig. 3a) 

Axial coordinate of lower limit of startline 
(ft or m, see Fig. 3a). (If ICON(3) point 
spacing option = 2 this value is recalculated 
using CORLIP(2)). 

Mach number (ICON(2)=0) or A/A* (ICON(2) = l) 
for startime 

Entropy of startline (cal/gm/°K or Btu/lbm/'^R) 
Mixture ratio (O/F) of startline 


Card 19a is used to input the gas startline information when the gas chemical 
equilibrium, frozen or ideal gas option is utilized in the solution. If a gas phase 
transonic solution is desired then this card reads in data for the transonic solu- 
tion. 36 
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Card 19- (Continued) 


Card 

Column 

Parameter 

Description 

19b* 

1-10 

CORLIP(2) 

Axial coordinate of upper limit of start - 
line (ft or m, see Fig. 3a) 


11-20 

CORRIP(6) 

Axial coordinate of lower limit of start - 
line (ft or m, see Fig. 3a) , 


21-30 

CORLIP(4) 

Mach number (ICON (2) =0) or A/A* 
(ICON (2) = 1) for startline 


31-40 

P 

Pressure for startline (atm) 


41-50 

T 

Temperature for startline (°R or °K) 

19c 

1-10 

CORLIP(2) 

Mass flow of nozzle (Ibm/ 8E10.6 

sec,gm/sec). If set = 0.0 

then the average O/F ratio 

(CORLIP(8))will be used to 

calculate the nozzle mass 

flow. 


11-20 

CORLIP(6) 

Initial station to start tran- 
sonic solution (ft, m) con- 
sistent with boundary 
equations. This station 
should be just upstream of 
the point in the nozzle inlet 
where it begins to contract. 


21-30 

' CORRIP(4) 

Axial position of nozzle 
throat (ft, m) 


31-40 

CORLIP(5) 

Axial station to stop the 
transonic solution (ft, m). 
Usually a station corre- 
sponding to an area ratio 
of around 3. 


41-50 

CORLIP(8) 

Engine mass average O/F 
ratio 

19d 

1-5 

lOFS 

Number of tables of O/F ratio 15 
versus radial position on the 
initial station (CORLIP(6)) 


(Maximum of 20). If IOFS=:l 
then constant O/F case -will 
be run using CORLIP(8) and 
card 19e is not used. 


This card is used to input the gas startline information when the gas chemical 
non- equilibrium option is utilized in the solution. 
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NOTE; Use card 19e when IOFS>l. Points should be input starting at nozzle 
centerline and going toward wall. This card inputs O/F distribution 
on the transonic initial station (COR LIP (6)). 


19e 

1-10 

OFST(I, 1) 

Radial position on transonic 8E10.5 

start line 


11-20 

OFST(I, 2) 

O/F ratio at above radial 
coordinate on start line 


21-30 

OFST(I+l,l) 

R 


31-40 

OFST(I+2,2) 

o/f 


41-50 

OFST(I+2,2) 

R 


51-60 

OFST(I+2,2) 

o/f 


61-70 

OFST (1+3,1) 

R 


71-80 

OFST(I+3,2) 

o/f 

NOTE; 

Use as many cards I9e as required to read in lOFS points. 

Card 20 

Startline Data 

F ormat 6E13 .7 

Do not use this 

card if lCON(2),:^2 or for gas-particle flow. Use feet 

if ICON(9)=0, meters 

ifICON(9) = l. Use Card 20a if ICON (1)<2. Use Card 20b 

if ICON(l)>2. 



Repeat this card in succession 

and in order of increasing R for 1=1, 2, . . . , 

ICON (3) 




Card 

Column 

Parameter 

Description 

20a 

1-13 

R 

Radial coordinate (R) of point I on 
startline (ft or m) 


14-26 

X 

Axial coordmate (X) of point I (ft or m) 


27-39 

EM 

Mach number at point I (dimensionless) 


40-52 

THETA 

Flow angle at point I (deg) 


53-65 

S 

Entropy at point I (cal/gm/°K or Btu/ 
Ibm/OR) 

- 

66-78 

OF 

Mixture ratio at point I (O/F) 

20b 

1-13 

, R 

Radial coordinate (R) of point I on start - 
line (ft or m) 

*This 

non- 

card is used to input the gas 
equilibrium option is utilized 

startline information when the gas chemical 
in the solution. 


Card 20a is used to input the gas startline information when the gas chemical 
non- equilibrium, frozen or ideal gas option is utilized in the solution. 
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Card 

Parameter 

Description 

14-26 

X 

Axial coordinate (X) of point L (ft or m) 

27-39 

EM 

Mach number at point I (dimensionless) 

40-52 

THETA 

Flow angle at point I (deg) 

53-65 

T 

Temperature at point I (°R or °K) 

66-78 

P 

Pressure at point I (atm) 


Card 21 

Cutoff Limits Data Format 8E10.6 

Required (See Fig. 3b) 

Column 

Parameter 

Description 

1-10 

CUTDAT(l) 

Radial coordinate defining upper limit of 
calculation regime (ft or m) 

11-20 

CUTDAT(2) 

Axial coordinate defining upstream cutoff 
limit (ft or m) 

21-30 

CUTDAT(3) 

Angle upper limit of calculation regime 
makes with horizontal (deg) 

31-40 

CUTDAT(4) 

Radial coordinate defining downstream 
cutoff limit (ft or m) 

41-50 

CUTDAT(5) 

Axial coordinate defining downstream 
cutoff limit (ft or m) 

51-60 

CUTDAT(6) 

Angle downstream cutoff line makes with 
horizontal (deg) 

Card 22 

Mesh Control 
Required (See 

Fromat 8E10.6 
Section 3.5.1 of Ref. 2) 

Column 

Parameter 

Description 

1-10 

STEP(3) 

Interior point insertion criteria (ft or m). 
See Section 3.5,1 of Ref. 2. 

11-20 

STEP(6) 

Axis point insertion criteria (ft or m). 
See Section 3.5.1 of Ref. 2. 

21-30 

■ STEP(9) 

Particle limiting streamline insertion 
criteria. 

31-40 

STEP(7) 

Point deletion criteria. 

41-50 

STEP(l) 

Prandtl-Meyer integration step size (deg). 

51-60 

STEP (8) 

Interpolation factor for calculating lower 
wall. 
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Fig. 3a ~ Startline Geometric Set-Up 



X, 


R, 

L 

X, 




© 


CUTDAT(l) 

CUTDAT(2) 

CUTDAT(3) 

CUTDAT{4) 

CUTDAT(5) 

CUTDAT{6) 


Normal terminated 
by downstream 
cutoffs . 

Normal terminated 
by plume boundary. 


NOTE: The normals must terminate on an upper 

boundary. Therefore, must have 

values such that the cutoff box will always 
be above the plume or solid boundary. The 
code will attempt to fill up the cutoff box 
with normals until fewer than six points 
remain on the normal. 


Fig. 3b - Cutoff Limits 
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Card 23 

Free -Molecular 
(This card is not 

Control Variables Format 6E10. 6 

used if ICON(13)=0) 

Column 

Parameter 

Description 

1-10 

VIBNO 

Reciprocal of the Knuds en number at which 
the vibrational energy mode thermally freezes 

11-20 

ROTNO 

Reciprocal of the Knudsen number at which 
the rotational energy mode thermally freezes. 

21-30 

TRANNO 

Reciprocal of the Knudsen number at which 
the translational energy mode thermally 
freezes. 

31-40 

CHARL. 

Characteristic length used m the Knudsen 
number calculation (nomally the nozzle 
exit radius) . 

41-50 

VISCC 

Reference viscosity (poise) if not input m 
thermo tables. 

51-60 

CON MM 

Viscosity relation temperature exponent 
if not input in thermo tables. 


Cards 24 through 36 are input only for two-phase solution. 


Card 24 

Column 

4 


5 


6-10 


Particle Solution Control Format 16I5 

(Use only if ISPECS >0) 


Parameter Value 
IZO 0 

1 

IWRITE 0 

(particle print 
flag) ^ 

2 

IDRAG 0 

1 


Description 

Nozzle wall equations are referenced to the 
nozzle throat. 

Nozzle wall equations are referenced to the 
nozzle exit plane. 

1 line of print for each particle (V, 9, AM, 
h, P, T) 

Above plus Re, AV, AT, viscosity, C^, Pr 

j^g, Nu/Nus, 

Use drag table coded in Kliegel program 
(Ref. 28). 

Use C. J. Crowe drag table coded internal 
to program (Ref. 29), 


All of above plus T , P , C_/C 
A -o ^ o’ o’ 
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Card 24 (Continued) 


Cojnmn 

14-15 


18 

24 


25, 30 


65 


Card 25 


Column 

1-10 


Card 26 


Column 

1-10 


Parameter Value 


Descript ion 


NSETS 


IPCHS 

JTEM(p 


JTEM(M) 

M = 2,ISPECS 


Particle 


0 Startline calculated by program 

(ICON (2)7^2) 

N Number of startline points at which 

particles are present for given startline. 
ICON (2) =2, N < ICON (3). Right adjust. 

0 No punch. 

1 Punch startline from transonic program. 

The elements of the JTEM(M)' arrary indi- 
cate which temperature/enthalpy table is 
to be used for particle species M. The 
value of JTEM(l) is always set equal to 
1 for particle species 1. 

0 Indicates that the particle species M 

temperature/enthalpy table will be the 
same as that for particle species 1. 

Cards 31, 32 and 33 are not required 
for particle species M. 

M Indicates that the particle species M 

temperature/enthalpy table will be input 
on Cards 31, 32 and 33 as Table M. 

N Indicates that the particle species N ' 
temperature/enthalpy table will be the 
same as that for particle species M. 
(N<M). Cards 31, 32 and 33 are not 
required for particle species N. 

Format 8E10.6 


Parameter 

XMASSP 


Description 

Ratio of particle total mass flow rate 
to gas mass flow rate. 


Particle Mass Flow Format 8E10.6 

Rate Fractions (Use only if 

ISPECS>0) 


Parameter Description 

PERTG(l) Ratio of particle No. 1 mass flow rate 

to total particle mass flow rate. 
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Card 26 (Continued) 

Column Parameter 

11-20 PERTG{2) 


Description 

Ratio of particle No. 2 mass flow rate 
to total particle mass flow rate. 


PERTG (ISPECS) 

Card 27 Particle Size Data 

(Use only if ISPECS >0) 

Column Parameter 

1-10 PSP(2, 1) 


Ratio of particle No. ISPECS mass flow 
rate to total particle mass flow rate. 

Format 8E10.6 


Description 

Radius of particle No. 1 (microns). 


Card 28 

Column 

1-10 


PSP(2, ISPECS) Radius of particle No. ISPECS (microns). 

Particle Mass Density Format 8E10.6 

(Use only if ISPECS >0) 

Parameter Description 

PSP(1, 1) Mass density of particle No. 1 (Ibm/ft^, 

or kg/m^) 


• • • 

PSP(1, ISPECS) Mass density of particle No. ISPECS 

(Ibm/ft^, or kg/m^). 

Card 29 Emissivity Data" Format 8E10.6 

(Use only if ISPECS >0) 

(e in Eq. (3.6) of Ref. 2) 

Column Parameter Description 

1-10 EMISS(l) Emissivity of particle No. 1. 


EMISS(ISPECS) 


Emissivity of particle No. ISPECS 


The emissivity and accommodation coefficients are used to determine the local 
energy exchange between the gas and particles via radiation. They normally 
produce negligible affects on solution and usually are set to 0 (zero). 
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Card 30 

Accommodation Coefficients Format 8E10.6 

(Use only if ISPECS>0) 

(a in Eq. (3.6) of Ref. 2) 

Column 

Parameter 

Description 

1-10 

ACC(l) 

Accommodation coefficient of particle 
No. 1. 

• 

ACC(ISPECS) 

Accommodation coefficient of particle 
No. 1 ISPECS. 

Card 31 

Particle Equation of State Format 4A6,I3, A6 

(Use only if ISPECS>0) 

Column 

Parameter 

Description 

1-24 

ALPHA 

Particle name (any name). 

28-33 

UNIT 

(Independent of 
ICON(9)) 

ENG Data input m English units 
MKS Use metric units 

Card 32 

Particle Data 

(Use only if ISPECS>0) 

Format 13, 12A6 

Column 

Parameter 

Description 

1-3 

NPTM(I) 

Number of temperature -enthalpy data 
points for this particle. If equal to 1, 
input liquid and solid heat capacities 
(see Card 33). Right adjust. 

Card 33 

Particle Enthalpy Data 
(Use only if ISPECS>0; 
units specified by Card 

Format 7E10.6 

31), 

Column 

Parameter 

Description 

1-10 

TM(I) 

Melting point temperature of particle 


No. I (°R in English units, °K in MKS 
unit s) . 


The emissivity and accommodation coefficients are used to determine the local 
energy exchange between the gas and particles via radiation. They normally 
produce negligible affects on solution and usually are set to 0 (zero). 
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Card 33 (Continued) 


Column 

11-20 

21-30 


Parameter 

HS(I) 

HM(I) 


Description 

Enthalpy of solid phase of particle No. I 
at melting point temperature (Btu/lbm 
or cal/gm) . 

Enthalpy of liquid phase of particle No. I 
at melting point temperature (Btu/lbm 
or cal/gm) 


If NPTM(I)=1, use following format. 


31-40 

41-50 

31-40 

41-51 

51-60 

61-70 


APHO(l, 1,1) Heat capacity of liquid phase of particle 

No, I (Btu/lbm-°R or cal/gm-°K). 

APHO(l, 2,1) Heat capacity ot solid phase of particle 

No. I (Btu/lbm -OR or cal/gm-°K). 

I'f NPTM(I) > 1 use following format. 

APHO(l, 1,1) Temperature for T-H table for particle 

No, I (OR or °K). 


APHO(l, 2,1) 
APHO(2, 1,1) 
APHO(2, 2, 1) 


Enthalpy for T-H table for particle No, I 
(Btu/lbm or cal/gm). 

Second temperature in T-H table for 
particle No. I (°R or °K) . 

Second enthalpy in T-H table for particle 
No. I. 


The above format (APHO(J, 1,1), APHO(J, 2,1)) is continued on successive 
cards of format 7E10.6 for J=l, 2, . . ,,NPTM(I). 

There are as many sets of cards 3l, 32 and 33 as there are different 
chemical species. 


Card 34 Input Startline Format 6E13.7 

(The following cards are required 
if ICON(2)=2 and ISPECS>0). 

Use Card 34a if ICON(l) < 2. Use Card 34b if ICON(l) > 2. 

Repeat this card lor 1 = 1, 2, . . . ,lCON(3) starting at point on nozzle axis. 


Column Parameter 

1-13 R 

14-26 X 


Description 

Radial coordinate of startline point I 
(ft or m) . 

Axial coordinate of startline point I 
(ft or m) . 


This card is used when gas chemical equilibrium, frozen or ideal gas option 
is selected. 
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Card 34a (Continued) 


Column 

Parameter 

Description 

27-39 

EM 

Mach number at point I. 

40-52 

THETA 

Flow angle at point I (deg) , 

53-65 

S 

Entropy at point I (Btu/lbm-°R or 
cal/gm-°K). 

66-78 

OF 

Gas total enthalpy (Btu/lbm or cal/gm). 

Card 34b* 

Column 

Paramet er 

Description 

1-13 

R 

Radial coordinate of startline point I 
(ft or m) . 

14-26 

X 

Axial coordinate of startime point I 
(ft or m) . 

27-39 

EM 

Mach number at point I. 

40-52 

THETA 

Flow angle at point I (deg) . 

53-65 

T 

Temperature at pomt I (°R or *^K) 

66-78' 

P . 

Pressure at point I (atm). 

Card 35 

Startline Particulate Data Format 15, 5X, 4E13.7 

(The following cards are required 
if ICON (2) =2 and ISPECS>0), 


Column 

Parameter 

Description 

5 

J1 

Particle number 

11-23 

HI 

Particle enthalpy at point I (Btu/lbm 
or cal/gm) . 

14-36 

RHOl 

Particle density at point I (slug/ft^ or 
kg/m^). 

37-49 

U1 

Particle axial velocity at pomt I (ft/sec 
or m/sec) . 

50-62 

VI 

Particle radial velocity at point I (ft/sec 


or m/sec) . 

Card 35 is repeated for each discrete particle size at each point on 
the start line where particles are present, starting at the nozzle wail and 
going toward the axis (reverse order of Card 34). 
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Card 36 Transonic Flow Data Format: Namelist 

(Use only if ISPEOO 
and ICON(2)/2) 

'Although, there are many parameters that may be input via the namelist 
DATA, most of these have already .been assigned values in the previous 35 
input cards; and some of the parameters do not apply to the transonic calcu- 
lation. Only those namelist parameters that could have a significant effect 
on the program are included below. The namelist data begins in Column 2 
with $DATA. The last card begins in Column 2 and contains only $END. 


Parameter 

Assumed Value 

THID 

Throat inlet half angle (deg) 

None 

THFD 

Fair mg angle (deg) 

(If THFD > THID no fairmg) 

5.0 

THJD 

Angle defining farthest downstream 
zone m transonic region (deg) 

9.0 

THIW 

Angle where start line intersects 
nozzle wall (deg) 

12.0 

RRT 

Throat wall radius of curvature 
divided by throat radius (>2.0) 

None 


ZAX Value of X where, startline None 

intersects nozzle axis, 
normalized by throat radius 
(If ZAX is not input the pro- 
gram will calculate a value.) 


(See Fig. 4 for an illustration of above parameters.) 


ZI 

Number of zones mto which the 
upstream portion of transonic 
zone is divided*. 

3.0 

ZJ 

Number of zones into which the 
downstream portion of transonic 
zone is divided* 

2.0 


47 


LOCKHEED • HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D7 84753 



Fig. 4 - Inlet eind Throat Parameters for a Gas -Particle Transonic 
Solution (see Ref. 28) 
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If ISKPY=7 read in cards 37, 38 otherwise leave out. 


Card 37 


Format: 15 

Format: 15 

Column 

Parameter 

Description 


1-5 

INUM 

Number of chemical species 
to be read in for this chemical 




system (maximum of 25) 



Card 38 

Card 38 reads in the species names for 
chemical system which does not appear 
in Table 1. 

18A4 

Column 

Parameter 

Description 


1-12 

(Left 

Adjusted) 

(SNAME(7,I, 
J), J-1,3) 

Species name (consistent 
with JANNAF 

3A4 

13-24 

I+l 


3A4 

25-36 

• 


3A4 

37-48 

• 


3A4 

49-60 

■ 


3A4 

61-72 

(SNAME(7, 
1+7, J), 

L 3) 


3A4 

Use 

as many cards 

38 as required to read in INUM species. 

If l^ISKPy^7 and ICON(l)=l card 39 is required, 

otherwise leave 

out. 




Card 39 

Species mole fraction card 
for SPF start line. 

Format: 8E10.6 

Column 

Parameter 

Description 


1-10 

SIFROZ(I) 

Species mole fraction 


11-21 

SIFROZ(I+l) 

Species mole fraction 



71-80 

SIFROZ(I+7) 

Species mole fraction 
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As many species mole fractions must be read in as there are species 
in the particular chemical system (see Table 1). If ISKPY=7, there should 
be INUM species. The order must be the same as the species appearing in 
Table 1 or were read in under card 38. Use as many cards 39 as necessary 
to read in all species mole fractions (25 maximum). 


Cards 40, 41, 42, 43 may be required only if a boundary layer calcu- 
lation is to be performed (IBL>0, card 5) otherwise they must be left out. 


Card 40 


Column 

1-5 


6-10 


11-15 


Boundary Layer Run Option Card Format: 

This card must be input if IBL>0 

parameter Value Description 

KR7 0 Turbulent boundary layer, 

calculation starts off 
laminar and will trip to 
turbulent when the Reynolds 
number based on momentum 
thickness exceeds 250. The 
program will select whether 
the flow is frozen or in equi- 
librium. 


515, 5X, 
E10.5) 


KR9 


1 Laminar boundary layer, 
frozen chemistry 

2 Laminar boundary layer, 
equilibrium chemistry 

3 Turbulent boundary layer 
if momentum thickness 
Reynolds number exceeds 
250, Equilibrium chem- 
istry. 

4 Same as 3 except frozen 
chemistry 

0 Adiabatic wall nozzle wall 

2 Specified wall temperature 

distribution. (Requires 
cards 43). 

4 Wall steady state energy 

balance 


ITPTS 


N where N is the number of 
wall points to input the 
nozzle wall temperature. 
This is used only if KR9=2 
(25 max.) 
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Card 40 (Continued) 


Column 

Parameter 

Value 

Description 


16-20 

ISPEC 

N 

N is the number of chemical 
species in the boundary layer 
edge gas. This is required 
only if ICON(l)=l (thermo 
data read from cards). 


21-25 

NUNIPP 

0 

Do not punch the BLIMPJ 
input data. 




1 

Punch the BLIMPJ input 
data. 


31-40 

XSTART 


Distance from the beginning 
of the boundary layer. Pro- 
gram defaults to l/2 a throat 
radius. 


Card 41 

Boundary Layer 

Edge Species Format: 

4(2A4, 

NOTE: 

This card is required only if ICON(i)=l. The order 
in which the species are read in does not matter. 

E10.6) 

Column 

parameter 


Description 


1-8 

SPECN(l) 


Chemical species name 
(Left Adjusted) 


11-20 

SPECMF(l) 


Species mole fraction 


21-28 

SPECN(2) 


2nd chemical species name 


31-40 

etc. 

SPECMF{2) 


2nd chemical species mole 
fraction 


As many card 41 are used as are 

required to read in ISPEC species, 

4 to a card. 




Card 42 

Transport Property Card 

Format: 

2E10.5 

NOTE: 

This card is required only if ICON(l) = l and the case being run is not 
a two phase case (ICON(9)=0). 

Column 

Parameter 


Description 


1-10 

XMU 


Gas viscosity in the -com- 
bustion chamber (poise) 


11-20 

PRER 


Turbulent Prandtl number 
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Card(s) 43 Nozzle Wall Temperature Distribution Format: 8E10.5 

NOTE: Card(s) 43 are only required if the wall temperature distribution option 
is selected (KR7=2). Input points starting from nozzle throat to the lip. 

Column Parameter Description 

I- 10 TWA{1, 1) Axial coordinate of wall 

location (ft, m) 

II- 21 TWA(1, 2) Wall temperature at above 

location (deg R, deg K) 

Input ITPTS pairs of location and temperature. Four pairs per card. 


52 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D7 84753 


Table 3 

MAGNETIC TAPE ASSIGNMENTS 
FOR TRAN72/RAMP/BLIMPJ PROGRAMS 


Unit Tape Unit Function 

1,8, 9j il> 12, 13 Internal Temporary Use Files 

2 BLIMPJ Output File 

3 Fiowfield Data Generated by RAMP 
Program (Output) 

4 TRAN72 Thermo Data (Input) 

10 Gas Properties Data Generated by 

TRAN72 program (Input) 
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2.3.2 Control Card Set-Up for UnLvac 1108 Exec 8 and Program Overlay 
Structures 

A typical run stream set-up for the Univac Exec 8 computer is pre- 
sented in this section. Also included is a table which gives the overlay 
structure for the RAMP program. 


NOTE: This schematic is typical of a run control scheme for the Univac 1108 
Exec 8 computer. It is presented to acquaint the user with magnetic 
tape and scratch area assignments. 


Table 4 gives the program overlay structure for the RAMP program. 
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flSG,T Rl T^F/Z/GO'^ 

UFF F'Rl T 

ASG^T niPtF///“0O 

uGr- TiBi P 

AFG-*T PPOO^F/Z/3^'1 

UFF A,PpOD 

AFGoT DlF! ,FZ/Z50P 

PC.C-- p,DU' 

AFG^T P.F 

ASG 1 T GAS 1 F 

USE IOiGAS 

ASG^T 1 1 »F 

ASG.T l?iF 

ASGiT PI I^^K///^PO 

USE 1 -3 » I N 

F!<FF TPPF 

A'=G»T T'-’FT. tF/Z/SPO 
ASG»TF STK,UnS»0/j 12V 
vOVE ST'< , 3 
COPY STK»PROo 
PE"IND STK 
C^PIN ST'<»TPF3> 

XUT A 

•K •)!*>. j!-,- DATA 
FH 
Fr<EF 

ASGiT TPFS.F///5P0 
PF."'rND STK 
POVE FTK«1 
COPIN STK,TPF'E 
CCPIM STK,TPFiS 
FI^EL STK 
XOT GO 

*****->t--^-'!-* k>A’^'‘P N0Z7LF 

XOT pUIMP 


. P-uN control card 


CORF AFSIGNMFNT 
cope- ASS T GMMFNT 

pr-'ocpAA" and pr-^oDucT data tape 

MOVE TAPE. TO 4 TH FILE 

COPY PRODUCTS DATA ONTO UNIT 4 

Rtv'IN’D PROGRAM TAPE 

COPY TRAN 72 PROGi’iAR INTO CORE 

execute TRAN 72 PROGRAM 

ERASE CORF 
CORF ASSIGNMENT 
CORF A SS' I GNMENT 
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■n 2 
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o ^ 

> 

pV P> 

10 "o 

c > 

> o 

C fTi 


XOT GO 

***« -K- 
FIN 
FIN 


No ULIMPJ J OAT 
1 Card of RA ip da 

R=^5T OF DATA WA 


• MOVE program tape TO SECOND FILE 
. COOY RAIMP CODE INTO COPE 
. COPY BLIMPJ CODE INTO COPE 
. RELEASE PROGRAM TAPE 

, execute ramp code FOP NOZZLE SOLUTION 
QATA -M- 

. execute BLIMPJ code for nozzle boundary I-AYFR 
A IT CAMF FROM RAMP OUTPUT FILE •«••!(•**■«• 

j EXECUTE RAMP E'CR PARTICLE DL CALC AND STARTLIne 

JA FOP PARTICL-E trajectory and exit startline 
s loaded on UNIT 13 IN FIRST RAMP EXECUTION **■«• 
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Table 4 

RAMP OVERLAY STRUCTURE 
Main 
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Main 




MAIN 




DRIVER 

REPEAT 


TRANS 

BLKDAT 

ALGINT 


PAR TIL 

BOUND 

VOFEM 


ONED 

EMOFV 

RITE 


JAMES 

IT SUB 

KIKOFF 


PROP 

POFEM 

UOFV 


FIND 1 1 

RHOFEM 

TKEY 


WDGI 

TEMTAB 

TOFV 


DCALC 

TOFEM 

IDMPFP 


LEGS 

PAGE 




ERRORS 

DRAGMR 



DROGCP 

SPCTX 



VOFEM 

TAB 



THERMO 

FABLE 



THERMl 

TOFH 



LAGRNG 

XSI 



SITER 




B1 




PLUMIN 




RGVOFM 

INITP 

IBOUND 


IDTAPE 

GASRD 

STARTV 


IDMXSI 

GASTAP 

WOFA 


IDMTAB 

IMPUT 

AOASTR 


PARTIN 

LIPIN 

SETHTG 


WTT 

WXANDR 

PART PH 


WG 

PAGVOF 

MASCON 


OUTPUT 

AXIS 

PLMOUT 


POINT 

LOGIC 



THERMT 

MAXT 



MASS 

PHI 



OFSET 

IDMPHI 
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SITER 

FCALC 

NEWT 

ABCALC 

CCALC 

PCALC 

STRMLN 
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B2 





PHASE 1 





PPATPT 

DOTPRD 




GAPPBI 

VEMAG 




SLDP 

SOKINT 




CHEM 

ROTERM 




EMOFP 

MAXTIM 




RGMOFP 

FNEWTN 




COEFEQ 

VMODl 




NEWENT 

VMODEL 




COEFF3 

INRSCT 




CHECK 







C3 



STRNOR 


OUT 



SOKSOL 


SDSKIP 

OUT BIN 


PHYSOL 


ENTROP 

VMOD2 




DELTAF 

POJH 


D1 

D2 

WEAK 

NORSCK 


CARCTR 

AVERAG 

ESHOCK 

TOFENH 


SEPLIN 

STGMOD 

TURN 





C4 



MCCRMK 


EXPCOR 



DECODE 


PAFH 

HYPER 


CODEH 


MOCSOL 

PRANDT 


CODEE 


OVEREX 

THETPM 


CODEF 





REDIST 





CZ 


C5 


C6 

MASSCK 


FREEMC 


IRAD 

THRUST 


PRFRBD 

WTFLOF 

ITERM 

INTEGR 


BOUNDA 

ITARM 

LIMITS 
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B4 


B5 

PUNEX 


BLMPIN 

ARASSL 


IDMTAB 

SPECIE 


SPECIB 

IDMPOP 


WALPRP 

POP 


TOB 

B6 



WTFLOP 



IDMPRO 



PRO 



IDMPDT 



PDT 



C7 

C8 

El 

PARLOK 

TRACEP 

DUMSYS 

INRSCP 

GAS 


TAPMOV 

START 


READF 



INTERP 
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2.3.3 Description of Unformatted Binary Output 

The binary tape output on unit 13 is described in this section. Initial 
input data are written on the first part of the data tape and gaseous and 
particle data are written out for each data point in the flow field. This tape 
is formatted so that it may be used by other auxiliary routines (plot, plume 
impingement or radiation). 

GROUP I — General Information 
Number of Records = 1 

Write { ) (HEADER(I),I=1, 60), ISPECS,IMETRIC 

• HEADER 

o run identification (2A4) 

• date (3A4) 

« description (55A4) 

• ISPECS = number of particle species to be 

considered 

• IMETRIC = 0 English flowfield units 

= 1 Metric flowfield units 


GROUP II - Gas Data 
Number of Records = 1 + IOF*IS 


Write { ) (BETA{I), 1 = 1,6), lOF, IS 

e BETA is gas identification name (6A4) 

• lOF number of total enthalpy cuts through 

"Mollier chart" (max = 10) 

« IS number of entropy cuts (max =2) 

DO M=l,IOF 
DO 1 = 1, IS 

Write ( ) IV,IDATA, ((TEMP(J,K), K = l, IDATA), J=l, IV),IVT, 
((CPM(J, K), K = l,3), J=1.IVT), R STAR, PIN F, EMIN F, 
GAMINF, FINF, EXINF, XSHIFT 

^Determined from ICON (9). 
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• IV 

• IDATA 

• GAMINF 

• IVT 

• RSTAR 
9 PINF 

a EMINF 
» EXINF 

• FINF 

• TEMP 

9 XSHIFT 


number of velocity cuts through ’’Mollier 
chart" for this total enthalpy and entropy 
+ 2 (max = 15) 

number of gaseous species present for this 
total enthalpy and entropy (max = 98) 

freestream isentropic exponent 

= IV-2 

throat radius (ft or meters) 

2 

ambient pressure (psf or Newtons/m ) 
freestream or external stream Mach number 
limit to which equation applies 

linear static pressure gradient (slope) 0 approach 

contains the following information for each value 
of lOF, IS 

nozzle length (ft or meters) 


12345678 _9 

1 p£ 


(;:) 


Htg P T S ijj 7 Me X, 

■ 


c 

c 

o 

o 

n 

n 

s 

s 

t 

t 

a 

a 

n 

n 

t 

t 


IVHtg P T S iji 7 M 


X, 


IDATA 


« • • « 


(iii 

Xp 


Mole 

Fractions 


( \ Species 
) Name 


• Pj 

• Htg 

A P 
a T 

a S 

a 


freeze pressure (atm) 

total enthalpy of the gas (cal/gm) 

pressure (atm) 

temperature (°K) 

entropy (cal/gm-°K) 

molecular weight (gm/gm-mole) 
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• y isentropic exponent 

• M chamber Mach number = 0 

c 

• M^' throat Mach number = 1 

• M Mach number for this table entry 

• CPM contains the following information 

12 3 

1 Pr A C 

2 P 

3 

IV T 

• Pr Prandtl number 

• H viscosity (poise) 

• C specific heat at constant pressure (cal/gm-°K) 

P 

GROUP III - Gas Particle Data 

Number of Records = ISPECS+1 

Write ( ) IDUM, ({PSP(I, J), 1=1,2), J = l, ISPECS) 

• IDUM dummy word 

,1 3 3 

• PSP(1,J) mass density of j particle (slug/ft or kg^/m ) 

« PSP(2, J) radius (ft or m) 

DO 1=1, ISPECS 

Write ( ) NTABl, TMELT, HSOL, HLIQ, (HFIT(N, 1, 1), HFIT(N, 2, 1), 

N=l, NTAB) 

« NTAB 1 number of table entries for this species 

• TMELT melt temperature (°R or °K) 

9 HSOL enthalpy before phase change (ft^/sec^-°R or m^/?sec^-°K) 

• HLIQ enthalpy after phase change (ft /sec - R or m /sec - K) 

9 HFIT(N, 1,1) temperature (°R or °K) 

• HFIT(N, 2,1) enthalpy (ft^/sec^ orm^/sec^) 

• NTAB number of table entries for this species 

Note that if NTAB =1 species is ideal and HFIT(1, 1, I) = CpL ( specific heat of 
liquid) and HFIT(1, 2, I) = Cp^ (specific heat of the solid). 
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GROUP IV - Flowfield Data 


Number of Records = 1 + 2 * ILAST 


Write ( ) (ILAST, 1=1,7), THRUST, AEXIT,IEXIT 

m ILAST number of data points on the following normal 

surface. If ILAST = 0 there is no information 
to follow 


e THRUST 
« AEXIT 

« lEXIT 


thrust (Ibj, or Newtons) 

^ 2 2 
exit plane area (ft or m ) 


exit flag 


j 0 if no exit 
1 1 if exit 


Write ( ) ((ITYPE, R, X, M, 9, S, /li , 6, Htg), 1=1, ILAST), (V, 1=1, ILAST), 
(Wg, 1=1, -ILAST), {(p,P,T,y,:^), 1=1, ILAST) 

• ITYPE identifies type of point (wall, shock, 

interior, etc.) 

0 input point 

1 interior point 

2 wall point 

3 free boundary 

4 upstream shock point 

5 Prandtl-Meyer point 

6 downstream shock point 

7 slip line 

8 shock intersection point 

9 vibrational mode frozen 

10 rotational mode frozen 

11 translational mode frozen 


e R 

• X 
« M 

• 0 

• S 

• M 

• 5 

• Htg 


radial coordinate (ft or m) 
axial coordinate (ft or m) 

Mach number 
flow angle (rad) 

entropy (ft^/sec^-°R or m^/sec^-°K) 
Mach angle (rad) 
shock angle (rad) 

gas total enthalpy (ft^/sec^ or m^/sec^) 
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• V 

• W 

g 

• P 

• P 

e T 

« y 

« H 

DO I = 1, IDAST 


velocity (ft/sec or m/sec) 

mass flow between this streamline and axis 
(slug/sec or kg^/sec) 

gas density- (slug/ft or kg /m ) 

z 

pressure (Ibf-ft or N/m ) 
temperature (°R or °K) 
isentropic exponent 

universal gas constant divided by molecular 
weight (ft^sec^-^R or m^/sec^-°K) 


Write ( ) (SPECN,I=1,NS) 

• NS number of gas species 

• SPECN species mole fractions 

Write ( ) (PO,I=l,ILAST) 

• PO pitot total pressure 


Output only for 
finite rate cases 


Output only if 
3 line printout 
selected ICON(8) 


DO I = 1, ILAST 

Write ( )ISP, ((U, V,T,H,p), J=1,ISP), ILIMIT 

number of particle sizes at this point 
axial velocity component (ft/sec or m/sec) 
radial velocity component (ft/sec or m/sec) 
temperature or °K) 


• ISP 

« U 

• V 

e T 

• H 

• P 

• ILIMIT 


2 2 2 2 
enthalpy (ft /sec orm /sec ) 

3 3 

particle density (slug/ft or kg^/m ) 

0| not a limiting streamline 

1 I is a limiting streamline 


NOTE; The flowfield data are repetitively stored on tape as indicated above - 
normal surface after normal surface. When ILAST = 0 the end of the 
data has been reached. 
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3. CONCLUSIONS AND RECOMMENDATIONS 


This document has described the new capabilities which have been 
added to the RAMP code. In addition to the basic capabilities of the RAMP 
code, which allow it to perform detailed two-phase flow calculations for 
sujf)ersonic nozzle/ plume flows, the following improvements have been made; 

® A nozzle exit plane start line can be punched or put on tape which 
will directly interface with the JANNAF Standard Plume Flowfield 
Code (SPF). This start line can include any or all of the following 
phenomena; two -phase flow effects, ideal, equilibrium or finite 
rate chemistry, variable O/F distribution, nozzle wall boundary 
layer, shockwaves and transonic effects, 

• A nozzle boundary layer solution can be performed with little or 
no additional input requirements. 

• The effect of the nozzle boundary layer on particles which enter 
the boundary layer can be considered. 

e A finite difference shock capturing numerical operator has been 
added. This option has been checked out for single phase ideal 
gas solutions only. 

9 A variable O/F transonic solution module has been incorporated 
into the code. This module will treat the transonic portion of 
nozzles for ideal and equilibrium chemistry for constant as well 
as variable O/f flows. 

e An improved two -phase transonic solution has been added which 
allows the treatment of small throat radii of curvature. At the 
present time this transonic calculations domain does not extend 
far enough downstream to generate a normal start line so that 
the throat radius of curvature must be greater than 1.5. This 
deficiency will be corrected and code modifications will be 
forwarded to the users. 

0 The code will generate an exit plane start line which includes 
the nozzle wall boundary layer. The subsonic portion of the 
boundary layer is expanded to Mach = 1.01 via the Seubold 
method so that a plume solution is possible. 

The updated version of the RAMP code has considerably more capa- 
bility than the original version. It will allow users to generate an exit 
plane viscous start line with little effort unlike the numerous hand calcu- 
lations and computer runs which were necessary in the past. 
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Additional work in the area of rarefied flows is necessary to make 
the plume calculations more applicable in the back flow regions of the 
plume. 


In its present form, the code has the capability of producing data for 
the following applications: 


e 

e 

s 

9 

e 

e 

e 


Gas/Gas -Particle Impingement (Heat-Transfer-Loads) 

Rocket Nozzle performance (Thrust, I ) 

s p 

IR Signatures (Radiating Species) 

RF Attenuation (Electron Densities) 

plume Radiation (Radiative Heat Transfer Gas/Particles) 
Vehicle Base Pressure 

Base Heating (Convection-Recirculation) 
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Appendix A 

EMPIRICAL INPUT DATA AND INPUT DATA SUGGESTIONS 
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Appendix A 

The results obtained from the Reacting and Multiphase Computer Pro- 
gram (RAMP) are very sensitive to data which are input. For two-phase 
cases the mean particle size, distribution of sizes, specific heats, mass 
density, particle melting temperature, chemistry assumptions and boundary 
equations are the primary input variables which determine the results. Each 
of these variables will be discussed in some detail and suggestions will be 
made as to what values can be used for aluminized propellants. It should be 
noted, however, that the data presented is not necessarily the best available 

Mean Particle Size 


Several different methods have been employed for obtainmg mean par- 
ticle size. Included are techniques which correlate mean size to throat diam- 
eter (Ref. A- 1), mean motor L (Ref. A-2) (chamber volume/throat area), 
chamber pressure, residence time, particle loading, maximum stable droplet 
size as well as combinations of each of these parameters. As a simple esti- 
mate of mean particle size the correlation of Delaney (Ref. A-1) based on throat 
diameter can be used: 


D = 4 d;^ (A.l) 

m t 

where D is the mean particle diameter in microns and D, is the throat 
m t 

diameter in inches. 


Particle Size Distribution 


For nozzle calculations in which no particle impingement on the wall is 
anticipated, one particle size at the mean size can be used. However, for 
plume calculations a knowledge of the particle size distribution is necessary. 
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Delaney in Ref. A-1 showed that the distribution of particles for smaller 
motors (D^ < 3.5 in.), followed a log normal distribution (Fig. A-1). For the 
large motors (D^ >3.5 in.) the data indicate that the size distribution follows 
a normal distribution (Fig. A-2), To use these distributions' move the curves 
up or down to the mean size at the 50% coordinate, then divide the curve into 
5 or 6 sections and determine the mean size that goes with each of these 
sections. Table A-1 gives an example of the size distribution (for 6 discrete 
sizes) which was determined from the curve in Fig. C-1. 

Particle Specific Heats, Enthalpies and Melting Temperature 

The values for particle specific heats, enthalpies and melting tempera- 
ture which the authors use are shown in Table A-E. The specific heats shown 
are used for the ideal approximation of particle enthalpy verses temperature 
(i.e., the specific heat for liquid and solid phases of the aluminum oxide are 
constant). The user may find the tables of temperature verses enthalpy in 
the JANAF Thermochemical Tables (Ref. A-3). 

Particle Mass Density 

The mass density for aluminum oxide is different for the solid and liquid 

phases.' Reference A-4 shows the mass density of liquid aluminum oxide (Af-O ) 

3 3 ^ ^ 

to bo 188 Ibm/ft . The mass density of solid Ai^^^ 250 Ibm/ft . For cases 

where the particle temperatures will be higher than the melting temperature 

for most of the flow field, the liquid mass density should be used. In cases 

where the particle temperature will be below the melting temperature (i.e , 

plumes) the solid mass density should be used. 

Chemistry Assumptions 

There are numerous chemistry assumptions which can be employed by 
the RAMP code. The various assumptions are: (1) ideal gas (constant specific 
heat ratio and molecular weight); (2) equilibrium; (3) frozen (constant molecular 
weight, varying specific heat ratio); (4) equilibrium/frozen (equilibrium with the 
molecular weight constant below a specified pressure) or (5) finite-rate chemistry. 
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Table A- 1 

LOG NORMAL PARTICLE SIZE DISTRIBUTION 
FOR HI 5 PC MOTOR 


Particle 

Diameter 

Percent Total 
Particle Mass 
Flow 

1.2 

10 

1.9 

20 

2.65 

20 

3,5 

20^ 

5.0 

20 

8.0 

10 


Table A-2 

THERMODYNAMIC DATA 


Liquid A^^O^ Specific Heat (C^ ) — 

.34 Btu/lbm°R 

Solid Ai^O^ Specific Heat (C^ ) — 

.32 Btu/lbm°R 

Enthalpy of Solid Phase of Af^O^ “ 

1358.9 Btu/lbm 

at Melting Temperature 


Enthalpy of Liquid Phase of Ai^O^ — 

1858.7 Btu/lbm 

at Melting Temperature 


Melting Temperature — 

4188°R 
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The type of chemistry assumption is very case dependent and also 
depends on the use of the flow field. Table A- 3 presents various cases 

and applications, along with suggestions as to the type of chemistry 
assumptions to be used. 

Finite rate cases can generally be started at the nozzle throat assuming 
the species distribution is in chemical equilibrium since this is valid for most 
propellant systems. Tables A-4 through A -6 present some reaction mech- 
anisms which may be used for applicable propellant systems. These reaction 
mechanisms were obtained from data presented in Ref. A-5. 

Boundary Equations 

The boundary equations which are mput to the code should be smooth 
and not contain discontmuities in either the slopes or coordinates where no 
discontinuities are physically present. Fictitious discontinuities can result 
in undesirable mass flow errors showing up during a solution. In Fig. A-3 
a description of the boundary equations for the nozzle throat and free boundary 
are presented. 

More complex nozzle contours may be input with discrete points which 
define the wall as a function of radial position and flow angle versus axial 
position. 
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Table A- 3 

SUGGESTED CHEMISTRY ASSUMPTIONS 
FOR VARIOUS APPLICATIONS 



Use specific heat ratio which exists at lip. 
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Table A-4 


REACTION MECHANISM FOR SPACE SHUTTLE SRM PROPELLANT 


H 

+Oh 

4-M 1 

=H20 

+ M1 

22 

6.10-26 2.0 

0.0 

0 

+H 

+M1 

= OH 

+M1 

21 

2.00-32 0.0 

0.0 

0 

+0 

+M 1 

= 02 

+M1 

24 

3.60-30 1 .0 

—340 • 0 

H 

+H 


= H2 

+M2 

22 

2.60-30 1 *0 

0*0 

OH 

+ ri 


= H2 

+0 

14 

1 .40-14-1 *0 

~ / 0 0 u • J 1 

OH 

-f-0 


= H 

+02 

1 1 

4.00-11 0»0 

0.0 j 

OH 

+H2 


= H20 

+h 

1 4 

1 .00-17-2.0 

-2900.0 

OH 

+OH 


= H20 

+0 

13 

1.00-1 1 0.0 

- 1 1 00.0 

CO 

+0 

+M3 

= C02 

+M3 

23 

2.00-33 0.0 

-4000.0 

OH 

+CO 


= C02 

+h 

14 

1. 10-19-2.0 

1600.0 

CL 


+M1 

= CL2 

+M1 

24 

4.30-31 1*0 

1 260 . 0 

H 

+ll 

+.VI4 

= hC.l 

+M4 

22 

3.00-30 1.0 

0.0 

OH 

+HCL 


= H20 

+CL 

1 4 

1 .00-14-1 .0 

-1000.0 

0 

+HCL 


= 0H 

+CL 

13 

2.00-12 0.0 

-4600 . 0 

H 



= HCL 

+CL 

13 

2.00-10 0.0 

-2400.0 

CL 

+H2 


= HCL 

+H 

13 

a. 00-1 1 0.0 

— 6 ^ 60 . 0 

-1 


Catalytic Species 


M^ = 3 H^O, 2 CO^; All others: 1.0 

M^ = 20 H, 10 H^O, 3 CO^, 2.5 All others: 1.0 

M 3 = 20 O^, 10 H 2 O, 3 CO^, 1.5 CO; All others: 1.0 

= 10 H 2 O, 5 HCi, 5 Ci, 5 H, 3 Cl^, 3 CO 2 ; All others: 1.0 
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Table A-5 


> 

t 


REACTION MECHANISM FOR ^2~^Z PROPELLANT SYSTEM 


ft 

+oh 

•fMl 

= H20 

4-M 1 

22 

6.10-26 2.0 

0.0 

0 

+H 

■H'U 

= OH 

■♦•Ml 

21 

2.00-.J2 0.0 

0*0 

o 

•fO 

•fMl 

=02 

■fMl 

24 

3 .60—30 1 e 0 

-340.0 

H 

+h 

-Ki'^2 

= H2 

-♦■M2 

22 

2.BO-30 1.0 

0.0 

OH> 

+H 


= H2 

-EO 

14 

1.40-14-1 .0 

- 7000.0 

OH 

-1-0 


= H 

■♦-02 

1 1 

O 

. 

o 

\ 

o 

o 

. 

0.0 

OH 

■+-H2 


= H20 

■fH 

14 

1 .00-17-2*0 

-2900.0 

OH 

-^0H 


=H20 

-fO 

13 

1.00-11 0.0 

- 1 100.0 


Catalytic Species 

Mj = 3 H^O; All others; 1.0 

= 20 H, 10 H^O, 2.6 All others: 1.0 
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Table A -6 

REACTION MECHANISM FOR LOX-RPi PROPELLANT SYSTEM 


I 

1 — • 
o 


H 

+OH 

+M1 

=H20 

+M'l 

22 

6. 10-26 2.0 

0.0 

O 

+H 

+M1 

= 0H 

+M1 

21 

2.00-32 0.0 

0.0 

O 

+0 

•HMl 

=02 

+M1 

24 

3.80-30 1.0 

-340.0 

H 

+H 

+M2 

= H2 

+M2 

22 

2.00-30 1.0 

0.0 

OH 

+ H 


= H2 

+0 

14 

1 .40-14-1 .0 

-7000.0 

OH 

+0 


= H 

+02 

1 1 

o 

. 

0 

1 

o 

o 

. 

0.0 

OH 

+H2 


= H20 

+H 

14 

1 .00-17-2.0 

-2900.0 

OH 

+0H 


= H20 

+0 

13 

1. 00-11 0.0 

-1 100.0 

CO 

+0 

+M3 

= C0<£ 

+M3 

23 

2.00-33 0.0 

-4000.0 

OH 

+Co 


= C02 

+ H 

1 4 

1 . 1 0-19-2.0 

1600.0 




Catalytic Species 

M j 

= 3 H^O 

2 CO 2 ' others: 1.0 

^2 

= 20 H, 

10 H 2 O, 3 CO^, 2.5 All others: 1.0 

M3 

= 20 O 3 , 

10 H 2 O, 3 CO 2 , 1.5 CO; All others: 1.0 
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RC != radius of curvature of the 
circular arc of the throat 
RT = throat radius 
XO - axial distance from the origin 
of the coordinate sjrstem to the 
thr oat 

0 ~ throat divergence angle corre- 
sponding to the maximum value 
for which the throat conic equa- 
tion applies 


The conic equation for this case would have the following form: 

A = -1 for an upper equation, -ti for a lower equation (-1 for this case) 
B = RC^ - XO^ 

C = 2X0 
D = -1 

E = -(RC + RT) 

Xmax = RC sinS + XO 


An example of a free boundary is shown in the sketch below. 



The free stream approach flow is 
inclined at 15 deg to the plume with 
a gamma [7) of 1.4, a Mach number 
of 10, and a static pressure of 0.1 
psfa. 


PIInF 

E 

gammainf 

MINF 

THETAINF 


O.i (psfa) 

0 (No pressure variation with axial distance) 

1.4 

10 

-15° 


Fig,A“3 - Sample of Boundary Equations 
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